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Abstract. Using the combined spectral and spatial resolving power of the Low Energy Transmission Grating (LETGS) on 
board Chandra, we obtain separate spectra from the bright central source of NGC 1068 (Primary region), and from a fainter 
bright spot 4" to the NE (Secondary region). Both spectra are dominated by discrete line emission from H- and He-like ions 
of C through S, and from Ee L-shell ions, but also include narrow radiative recombination continua (RRC), indicating that 
most of the observed soft X-ray e mission arises in low-temperature (kTe ~ few eV) photoionized plasma. We confirm the 
conclusions of Kinkhabwala et al. (2002b), based on XMM-Newton Reflection Grating Spectrometer (RGS) observations, that 


the entire nuclear spectrum can be explained by recombination/radiative cascade following photoionization, and radiative decay 
following photoexcitation, with no evidence for the presence of hot, collisionally ionized plasma. In addition, we show that this 
same model also provides an excellent fit to the spectrum of the Secondary region, albeit with radial column densities roughly 
a factor of three lower, as would be expected given its distance from the source of the ionizing continuum. The remarkable 
overlap and kinematical agreement of the optical and X-ray line emission, coupled with the need for a distribution of ionization 
parameter to explain the X-ray spectra, collectively imply the presence of a distribution of densities (over a few orders of 
magnitude) at each radius in the ionization cone. Relative abundances of all elements are consistent with Solar abundance, 
except for N, which is 2-3 times Solar. Finally, the long wavelength spectrum beyond 30 A is rich of L-shell transitions of Mg, 
Si, S, and Ar, and M-shell transitions of Fe. The velocity dispersion decreases with increasing ionization parameter, which has 
been deduced from the measured line intensities of particularly these long wavelength lines in conjunction with the Fe-L shell 
lines. 
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1. Introduction 

The bright Seyfert 2 galaxy, NGC 1068, has been studied ex¬ 
tensively for many years at optical, UV, IR, radio, and X-ray 
wavelengths. In the unified mode l of active galactic nuclei 
(AGN) (Miller & Antonucci 


1983 


Antonucci & Miller 


1985 
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Recently, Kinkhabwala et al. (2002b) presented the X-ray 


1993), the soft X-ray spectra of Seyfert 2 galaxies 


Antonucci 

are expected to be strongly affected by emission and scattering 
in a medium irradiated by the nuclear continuum. However, 
soft X-ray emission can also be produced by collisionally- 
heate d gas associated with shocks in starburst regions (Wilson 
et al. 1992). High resolution X-ray spectroscopy, now becom¬ 


ing available with the grating spectrometers on Chandra and 
XMM-Newton, can provide a means of distinguishing between 
these two interpretations. 


spectrum of the nuclear regions in NGC 1068 obtained with 
the Reflection Grating Spectrometer (RGS) on XMM-Newton. 
They show that the soft X-ray emission is dominated by lines 
from H-like and He-like ions of all low-Z elements from C 
through Si, as well as a complex of L-shell lines from Ne- 
like through Li-like Fe. Most of the lines are significantly 
blueshifted and broadened, with characteristic velocities ~ sev¬ 
eral hundred km s“^. Of particular interest is the presence of 
strong and narrow RRC of C, N, O, and Ne ions, which im¬ 
ply that most of the observed soft X-ray flux arises in low- 
temperature {kT few eV) plasma. There is also excess emis¬ 
sion (relative to pure recombination) in all higher series reso¬ 
nance lines of the H-like and He-like species, which is shown 
to be a result of direct photoexcitation by the continuum. 


Kinkhabwala et al. ( |2002b[ ) further present a new quanti¬ 
tative model for the nuclear regions of NGC 1068 that pro¬ 
vides an excellent fit to their observed spectrum. The model 
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involves a cone of warm plasma photoionized and photoex- 
cited by the same nuclear power-law continuum (Kinkhabwala 
et al. 2002c). In this model, each ionic line series is ht inde¬ 


pendently, with only three free parameters; the radial column 
density for that ion, the velocity width, and the overall normal¬ 
ization (the product of the nuclear luminosity and the covering 
factor). Interestingly, the derived column densities are consis¬ 
tent with those measured using absorption features in high res¬ 
olution X-ray spectra of Seyfert 1 galaxies, thereby providing 
strong support for the unihed model. 

Given the 15" spatial resolution of the XMM-Newton tele¬ 
scopes, it was not possible with the RGS data to derive any 
useful information on the spatial dependence of the observed 
emission line spectrum. This issue is best addressed with 
the transmission grating spectrometers on the Chandra X-ray 
Observatory. In this paper, we present the data acquired with 
the LETGS (see als o the c ompanion paper presenting HETGS 
data by Ogle et al. 2002 ). Earlier Chan dra imaging observa¬ 
tions of NGC 1068 (Young et al. 2001 ) had shown the pres¬ 
ence of a bright compact source, ~ 1".5 or 108 pc in ex¬ 
tent (1"= 72 pc), taking the d istanc e to NGC 1068 to be 
14.4 Mpc (Bland-Hawthorn et al. 19971. This is coincident with 
the brightest radio and optical emission, together with an ex¬ 
tended cone of emission (out to 7"; 504 pc) pointing to the 
NE, which coincides with the NE radio lobe and gas in the nar¬ 
row line region. A discrete region of emission is observed in 
the cone at about 4" (288 pc) from the bright compact source. 
We show here that the spectrum of this region differs signih- 
cantly from that of the compact central source, but that it is also 
well described by the same radiation-driven model presented 
by Kinkhabwala et al. (2002b), albeit with lower radial column 
densities. The lower column densities are expected, and are a 
natural consequence of the larger distance of this extended gas 
from the source of the ionizing continuum. 

The LETGS also has a signihcantly larger bandpass than 
the RGS and allows us to investigate the nature of the spectrum 
down to the region of the Ee K-shell transitions near 2 A, and 
out beyond 35 A, where we detect L-shell transitions of Ar, S, 
Si, and Mg. 


2. Observation and Data Analysis 

Since the extended nature and the orientation of the X-ray 
emission from NGC 1068 were known prior to the time the 
grating observations were planned, the roll angle was selected 
such that the dispersion angle was nearly perpendicular to the 
NE extent, making it possible to separately study the emission 
from the central nuclear component (hereafter, the Primary re¬ 
gion), and a bright spot in the extended ionized cone (hereafter, 
the Secondary region). The HETGS (Ogle et al. 20021, and the 
LETGS (this work) observed the source on December 4 and 5, 
2000, for 46 ks and 77 ks, respectively. The LETGS conhgura- 
tion for this measurement involved the ACIS-S as the readout 
detector, in place of the usual HRC. We used CIAO version 
2.1 to generate the event hie. In order to verify the wavelength 
scale as generated with CIAO, we ran a Capella calibration ob¬ 
servation (OBSID 00055), and did a recursion analysis based 
on a few strong spectral lines. There appears to be a very small 



Fig. 1. LETGS zero order image. North is up and East is left. 
The strong bar (top left-bottom right) represents ’’out of time 
events” from the CCD-readout. The beginning of the spectral 
image, minus order bottom left and positive orders, top right 
can just be seen. Notice the dispersion direction is close to per¬ 
pendicular to the extended NE-direction. The point source in 
the lower right quadrant has a very hard spectrum. The angular 
size of the total image is 60" x 75". 


systematic offset proportional to the wavelength, which is close 
to a similar systematic offset found during the an alysis of the 
HETGS measurement of NGC 4151 (Ogle et al. |2000| ). The 
correction is -f 0.7 mA per A. 

The LETGS zero order image is displayed in Eig. The 
diagonal streak in this image is not the dispersed spectrum, it 
arises from “out of time” events in the CCD-readout. The short 
wavelength edges of the spectrum can be discerned as spots in 



Fig. 2. Zero-order X-ray intensity contours plotted over a n op¬ 
tical [O III] (5007 A) HST image (Bruhweiler et al. 2001 ). The 
Primary (“P”) and Secondary (“S”) extraction regions are indi¬ 
cated. The X-ray/optical correlation is remarkably good. 
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Fig. 3. Intensity as a function of the cross dispersion direction. 
The Primary (“P”) and Secondary (“S”) regions are indicated. 


the lower left and upper right of the figure. Note that the source 
extent is nearly perpendicular to the dispersion direction. In 
Fig. this X-ray image is shown as a contour plot overlayed 
on an [O III] optical image obtained by HST. As can be seen, 
the X-ray extension correlates well with the extended [Olll] 
emission. We have indicated in Fig. ^ the source extraction re¬ 
gions, which we use to isolate the spectrum of the Primary and 
Secondary regions. 

A spatial cut through this zero order image in the cross¬ 
dispersion direction is plotted in Fig. ||. Clearly visible on the 
left side (NE-direction) is the small Secondary peak, which 
corresponds to our Secondary region. Here again, the extrac¬ 
tion regions for the Primary and Secondary spectra are indi¬ 
cated. The boundary between them was chosen based on the 
O VII line ratio variation along the cross-dispersion direction 
(see Sect. 3). 

The zero order profile in the dispersion direction of the 
Primary region is well represented by a Gaussian with the 
centroid at —0.0056 ± 0.0007 A, and with a width given by 
cr = 0.0317 ± 0.0008 A. This profile is slightly broader than 
the instrument profile for a point source (FWHM = 0.06 A; 
a = 0.0255 A), and is therefore affected by the geometry of 
the source extent. Similarly, the profile of the Secondary re¬ 
gion in the dispersion direction is represented by a Gaussian, 
with the centroid at -1-0.032 ± 0.002 A, and with a width 

cr = 0.062 ±0.002 A. 

We checked for time variability in the zero order flux. As 
is to be expected for this spatially extended source, we find no 
significant variations on timescales in the range ~ few seconds 
up to the length of the observation. 

3. Spectral Analysis 

We present spectra of the Primary and Secondary regions in 
Figs. Q and H respectively. We combine m = — 1 and m = ±1 
orders to obtain these spectra. Both spectra are plotted over the 
full range of the LETGS of ~ 1.5 — 60 A. 


Table 1. He-like R ratios for Primary and Secondary region 
spectra. 


Ion 

Primary Secondary 

Cv 

N VI 

Ovii 

NeIX 

11 ±17 1.5 ±1.1 

3.3 ±1.5 — 

5.4 ±0.9 4.5 ±2.0 
2.3 ±0.5 1.4 ±0.6 


Table 2. He-like G ratios for Primary and Secondary region 
spectra. 


Ion 

Primary 

Secondary 

Cv 

2.5 ± 1.4 

1.8 ± 1.3 

N VI 

3.0 ± 1.1 

1.1 ±0.9 

Ovii 

4.6 ±0.9 

1.6 ±0.4 

NeIX 

2.0 ±0.3 

1.4 ±0.4 


The Primary region contains approximately ninety percent 
of the flux (”P” region in Pigs. ^ and |^). As is expected, this 
spectrum closely resembles the RGS spectrum (Kinkhabwala 
et al. 20021), which contains all of the soft-X-ray flux from 
NGC 1068. In addition to the H-like and He-like series of C, 
N, O, Ne, Mg, and Si and Pe L-shell transitions (from Pe XVII- 
Pe XXIV), which were observed by the RGS, the PETGS also 
extends down to Pe K at 1.8 A, and, at longer wavelengths, 
from the He-like C triplet at 40 — 42 A (which is just outside 
the RGS range) out to 60 A, resolving a forest of lines from 
P-shell transitions in mid-Z ions such as Ar, S, Si, and Mg. We 
also note that lines in the ^ 5 — 10 A region of the specrum, 
though detected by the RGS, are better resolved by the PETGS. 

As was also found in the RGS spectrum, the PETGS re¬ 
solves several prominent RRC, which, in addition to the rela¬ 
tive strength of the forbidden line in the He-like triplets, pro¬ 
vides the main evidence for recombination in low temperature, 
photoionized plasma. Many higher order transition lines, ob¬ 
served by the RGS, are also resolved by the PETG S; thes e lines 
are indicative of photoexcitation (e. g., Sako et al. 2000 ). 

Overall, many features observed in the Primary spectrum 
are also observed in the spectrum of the Secondary region (”S” 
region in Pigs. ^ and ||), which comprises approximately ten 
percent of the total flux from NGC 1068. 

One main difference between the two spectra is in the He- 
like triplet ratios, R = f /i and G = {i + f)/r, which are listed 
for several He-like ions in the Primary and Secondary spectra 
in Tables [| and 

The G ratios in the Primary spectrum are inconsistent with 
either pure recombination or collisional ionization equilibrium, 
suggesting the presence of photoexcitation , as wa s concluded 
for the RGS spectrum (Kinkhabwala et al. 2002b ). The G ra¬ 
tios in the Secondary region are consistent with hot collisional 
plasma or with photoexcitation in a recombining, photoionized 
plasma. 

The R ratios are overall less useful for distinguishing be¬ 
tween pure recombination and collisional ionization equilib¬ 
rium. They are also more difficult to measure due to the weak¬ 
ness and/or blending of the intercombination line in most 
triplets. However, we note that the well-measured R ratio for 
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g. 4. Spectrum of the Primary region between 1.5 and 60 A, in five panels to cover the full range. The bin size is equal 
025 A. Indicated are some of the strongest lines. All panels are the positive and negative orders taken together. 
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Fig. 5. Spectrum of the Secondary region, wavelength scales as in previous figure for easy comparison. 
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Fig. 6. Forbidden line to resonance line ratio of the O Vll triplet 
as a function of cross dispersion angle. The bins along the cross 
dispersion axis are one arcsec wide; negative numbers are in 
the NE direction, positive numbers in the SW direction and bin 
zero is centered at the maximum intensity. One clearly notices 
an abrupt change over from a value of 3.6±0.6 for the Primary 
region to 1.0±0.3 for the Secondary region, around the bound¬ 
ary between these regions at —2". 


O Vll is significantly greater than expected for p ure rec ombi- 
nation, as was also found in Kinkhabwala et al. ( 2002b ), who 
suggested that this may be due t o inner-shell ionization of O VI 
(e.g., Kinkhabwala et al. 2002a ). 

Since the O Vll triplet is the strongest feature in both spec¬ 
tra and since the line ratios are different for both spectra, we 
measured the forbidden-to-resonance line ratio //r as a func¬ 
tion of cross-dispersion, in steps of 1", see Fig. ^ There are 
clearly two distinct regimes, the Primary region with a high ra¬ 
tio of //r = 3.6 ± 0.6, and the Secondary region with a value 
of 1.0 ± 0.3. This changeover occurs within a spatial extent of 
less than one arcsec. 


3.1. Fits Using an Irradiated Cone Model 

We find that the H-like and He-like ionic line series observed 
in the Primary and Secondary region spectra are fit well with 
a simple irradiated cone model, as was us ed for t he RGS spec¬ 
trum of NGC 1068 (Kinkhabwala et al. 2002b ). The cone is 
irradiated at its tip by an inferred nuclear power-law contin¬ 
uum. Though our view to the nucleus is completely obscured 
by intervening material (dusty torus), our perpendicular view 
of the irradiated cone is unobscured, allowing for observation 
of recombination/radiative cascade following photoionization, 
and radiative decay following photoexcitation. This simple ir¬ 
radiated cone model has been incorporated into XSPEC as the 
local model photo (Kinkhabwala et al. 2002a), which has been 
used to generate the fits to the Primary and Secondary regions 
as presented below. 

We assume that the overall normalization, fLx, which is 
the covering factor / = Q./Att times the power-law luminos¬ 


ity Lx is the same for all ions. For ease of comparison of the 
LETGS results with the RGS results and between the Primary 
and Secondary region spectra, we assume the same value for 
the normalization of JLx = 10“^^ ergs s“^. 

We also assume a single radial velocity distribution spec¬ 
ified by Gaussian width The velocity width needed to 
explain the Primary region spectrum of = 100 km s“^, 
however, is a factor of two smaller than that required to explain 
the RGS spectra. For simplicity, we use this velocity width for 
the Secondary region as well. 

The only free parameters left are the individual radial ionic 
column densities. We allow for individual line absorption, but 
assume that photoelectric absorption is negligible. With this as¬ 
sumption, we can simplify the calculations by irradiating each 
ionic column density separately with an initially unabsorbed 
power law, since the effect of individual line absorption due 
to all other ions, which removes only a small fraction of the 
flux, can be safely neglected. Alternatively, using the method 
employed by photo, we can irradiate all ionic column densities 
together, assuming all ions have the same radial distribution. 
The radial ionic column densities we infer below are indeed 
consistent with the assumption of negligible photoelectric ab¬ 
sorption, and we have checked that either method yields the 
same result. Moreover, as we argue in Sect. there likely ex¬ 
ists a broad range of densities at each radius, implying that the 
assumption that all ions have a similar radial distribution may 
also be physically correct. 

We demonstrate below that these simple model assump¬ 
tions allow for a good fit to both the Primary and Secondary 
region spectra. 


3.1.1. Primary Region 


Using the above-quoted global values of JLx = 10^^ ergs s“^ 
and = 100 km s“^, we have fit the individual radial ionic 
column densities. The best fit values are given in Table |], and 
the final fit using these values is shown in Fig. ^ These column 
densities are very similar to the results obtained for the RGS 
spectrum. 


The ratio of the ground state RRC to the 2p^ 1 s transitions 
for each ion rules out the presence of any additional hot, col- 
lisional plasma. The forbidden lines in the He-like triplets of 
N, O, and Ne, however, appear significantly stronger than pre¬ 
dicted. This is likely due to inner shell ionizations in the cor¬ 
responding Fi-like ions, which enhance only the forbidden line 
in the triplet (for further explanation, see, e. g., Kinkhabwala 
et al. 2002a ). This observation and proposed mechani sm was 
also noted in the RGS paper (Kinkhabwala et al. 2002b). 


3.1.2. Secondary Region 

Using the same global values as above, we have fit the column 
densities for the Secondary region. The best fit values are given 
in Table |^, and the final fit is shown in Fig. |[ The column 
densities for the Secondary region are a factor of a few less 
than the column densities inferred for the Primary region. 
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Fig. 7. Fluxed spectrum of the first order of the Primary region, together with the model. The model does not yet include Fe 
L-shell transitions. 


The explanation for this enhancement of the resonance line 
for lower column density is simple (Kinkhabwala et al. 2002b). 
Keeping cr™'^ constant, at low column densities the resonance 
line is dominant in the He-like triplet, since this transition has 
a high oscillator strength and is unsaturated. At higher column 
densities this transition begins to saturate, and recombination 
begins to increase in significance (since the photoelectric edge 
saturates only at high column density). Generically, then, at low 
column densities, the resonance line is relatively strong, but at 
high column densities, it saturates, becoming relatively weaker 
compared to lines due to pure recombination. 


The ratio of the RRC to the 2p^ls transitions again rules 
out the presence of a significant additional hot, collisional 
plasma. Though, given the lower quality of the Secondary spec¬ 
trum compared to the Primary spectrum, this conclusion is 
correspondingly less robust. However, the presence of strong 
higher order lines (np—>ls), such as N VII Ly 7 and O VII HeJ, 
which can only be explained by photoexcitation, strengthens 
the argument against an additi onal ho t, collisional plasma con¬ 
tribution (Kinkhabwala et al. 2002b i. The formal upper limit 
to the emission measure of any hot, collisional plasma in the 


0.2-1 keV temperature range is 10®^ cm ^ for the Secondary 
region and 2 x 10®^ cm“^ for the Primary region. 

We also note that the temperatures for the RRC are the same 
as for the Primary region (in particular, compare the RRC of 
O VII in both spectra). 


3.2. Elemental Abundance Estimates 


In order to obtain an estimate of the relative elemental abun¬ 
dances, we derive a total hydrogen column density from the 
individual ionic column densities assuming Solar abundances. 
This hydrogen column density A^h is a function of the ion¬ 
ization parameter (0, according to the ^ in which each ion 
forms. For this analysis, we employ the same set of photoion¬ 
ization balance calculations from XSTAR (Kallman & Krolik 


Kallman & Krolik 1995) as used in the analysis of the Seyfert 1 
galaxy NGC 5548 (Kaastra et al. 2002) and assume that the ob¬ 
served flux for each ion originates predominantly from plasma 
at an ionization parameter where the relative concentration of 
that ion reaches a ma ximum . Adopting Solar abundances from 
Anders & Grevesse ( 1989 ), except for the recently improved 
measurement of the Solar O abundance taken from Allende 
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Fig. 8. Spectrum of the Secondary region, together with the model. 


Table 3. NGC 1068 Ionic column densities of Primary and 
Secondary regions. Italicized values for the temperatures were 
taken arbitrarily for ions which do not show clear RRC. 


Ion 

Nion Primary 

Nion Secondary 

kT (eV) 

Cv 

~TW^ 

3.10^' 

2.5 

Cvi 

8.10^^ 

4.10^'^ 

4 

Nvi 

3.10^^ 

1.10^'^ 

3 

Nvil 

5.10^^ 

LlQi" 

4 

Ovil 

7.10^^ 

2.10^'^ 

4 

0 VIII 

6.10^^ 

2.5.10^'^ 

6 

NeIX 

2.10^^ 

8.10^® 

10 

Nex 

TlQi" 

6.10^® 

10 

Mg XI 

TlQi" 

2.5.10^® 

10 

Mg XII 

TlQi" 

2.10^® 

10 

Si XIII 

1.1.10^^ 

3.10^® 

10 

Si XIV 

1.1.10^^ 

2.10^® 

10 

Sxv 

4.10^® 

1.10^® 

10 

Sxvi 

6.10^® 

1.10^® 

10 


Prieto et al. ( 2001 ), we obtain for the Primary compo¬ 

nent as plotted in Fig. 

Note that these column densities are consistent with the es¬ 
timated column of log TVh = ~21.9 derived from the L-shell 


lines (Sect. ||). As explained in detail in that section, this value 
of A^h is consistent with the observed line fluxes measured for 
the Fe lines assuming that they arise exclusively from photoex¬ 
citation. However, since most of the bright Fe lines are satu¬ 
rated, this estimate is very uncertain and thus of less value for 
abundance measurements. We estimate the accuracy of the col¬ 
umn densities based on the K-shell lines to be within ~ 25 % 
(relative to each other). The error bars in Fig. ^ reflect these 

The data point corresponding to H-like Ne in Fig. ^ seems 
to be somewhat low compared with the general trend suggested 
by the other ions. Therefore we consider this point on the plot 
less trustworthy for an abundance analysis. This inconsistency 
could be due to the Lya line of Ne X being contaminated by 
Fe XVII. 

Deviations from Solar abundances in Fig. ^ would appear 
as clear vertical offsets between elements that cover similar re¬ 
gions of However, it can be seen that obvious offsets between 
elements are not seen. In fact, the column density Nn does not 
vary by much in the observable range of This implies that 
for the most part the relative abundances (no measurement is 
available for H) are consistent with Solar, or at least that there 
are no conclusive deviations, of more than about ±50%, from 
Solar abundances. On the face of it, this is in stark contrast with 
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Fig. 9. Derived hydrogen column density as a function of ion¬ 
ization parameter for C, N, O, Ne, Mg, Si and S. 


previous moderate-resolution X-ray observations of NGC 1068 
as well as with UV and optical measurements from the nar¬ 
row line region. Netzer & Turner (1997) as well as Marshall 
et al. (1993) have fitted the ASCA and BBXRT X-ray spectra, 
respectively, of NGC 1068 with models that require an Fe/O 
abundance ratio of more than 8 times Solar. This could be at¬ 
tributed to the severely incomplete Fe-L model that was used 
in those works. Also, photoexcitation was not treated explicitly 
as shown to be necessary in Kinkhabwala et al. (2002b I. In the 
ASCA case, the weak signal in the low-energy region of the O 
features might also be to blame. On the other hand, the result of 
those papers that all other abundances except for Fe and O are 
Solar is consistent with the present measurement, with the ex¬ 
ception of N as discussed below. Nitrogen apparently was not 
included in those earlier models. 

According to Fig. ^ N is the only element that seems to lie 
consistently above the data points of other elements (namely 
O and C) with which it overlaps in ^ space. The N abundance 
enhancement found here is of the order of a factor 2 —3 com - 
pared to C and O (also observed in Kinkhabwala et al. 2002b ). 
This could be due to an enrichment of the AGN medium by 
WR type stars close to the nucleus. An enrichment of nitro¬ 
gen has also been found in other AGN, for example in IRAS 
13349H-2438 (Sako et al. 2001). Finally, we note that for the 
Secondary region, we find very similar results only that the A^h 
is about 10^^ cm“^, which is 3—5 times smaller than for the 
Primary region. Additionally, the error bars for the Secondary 
region are obviously larger than for the Primary region. 


4. Line Broadening and Line Shifts 

4 . 1. Line Broadening 

We have used a number of strong emission lines, the Lya lines 
of Si, Al, Mg, Ne, O, N and C, and the forbidden lines of the 
He-like triplets of O and N, to look for velocity broadening 
(Table Q. All lines of the Primary region appear to be broader 
than the zero order profile, although the formal errors are large. 
The most reliable velocity numbers are obtained from O viii 



Velocity (km/s) 


Fig. 10. Spectrum of the O VII forbidden line in -f 1 order for 
the Primary region. The velocity scale is set with respect to 
the rest frame of NGC 1068. The solid line shows the best fit 
model with an average blueshift of —250 km/s and a velocity 
broadening of 320 km/s. The dashed line is the best fit with no 
velocity shift and broadening. Note that in the latter case the 
line profile is given by the zeroth order profile of the Primary 
as discussed in Sect. 3.1. The bin size in wavelength space is 
25 mA. 



Velocity (km/s) 


Fig. 11. Spectrum of the O VII forbidden line in -|-1 order for 
the Secondary region. The velocity scale is set with respect to 
the rest frame of NGC 1068. The solid line shows the best fit 
model with an average redshift of 4-360 km/s and no velocity 
broadening. The dashed line is the best fit with no velocity shift 
and broadening. Note that in the latter case the line profile is 
given by the zeroth order profile of the Secondary as discussed 
in Sect. 3.1. The bin size in wavelength space is 37.5 mA. 


Ly a and the forbidden line of the O Vii triplet, cr^ being 600 
and 320 km/s, respectively. The Si, Al, and Mg Lya lines show 
values of about 1000 km/s but with large errors. Fig. shows 
the O VII forbidden line of the Primary region in velocity space. 

For the Secondary region the measured widths of the lines 
are consistent with zero, with upper limits of about 1000 
km s“^. 
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Table 4. Line strengths and velocity shifts of some of the stronger lines in the Primary and Secondary region. Line strengths are 
at the source and in units of 10^'^® photons s“^. We adopt a galactic absorption iVn = 2.95 x 10^° cm“^ (Murphy et al. 1996| ) 
and a redshift of z = 0.00379 ± 0.00001 (Huchra et al. 1999| ). The errors quoted are statistical errors; the systematic error in the 
velocities in columns 5 and 6 range from 100 to 30 km s“^ for Si and O, respectively. All velocity values are quoted with respect 
to the systemic velocity. The line strengths and Av values of the Ne X Lya line are probably contaminated with Fe-L lines and 
are therefore considered less reliable. 


Line ID 

Pred. A 

A 

Primary 

photons (lO'^^s”^) 

Secondary 
photons (lO'^^s”^) 

^Vprimary 

km 

^^Secondary 

km 

Si XIV Lya 

6.180 

0.22±0.03 

0.05±0.02 

- 940±290 

-1-470 ±1650 

Mg XII Lya 

8.419 

0.25±0.03 

0.12±0.03 

- 1000±360 

-1-290 ±1030 

Ne X Lya 

12.132 

0.83±0.11 

0.50±0.08 

- 550±150 

-1- 1450 ±590 

0 VIII Lya 

18.969 

2.95±0.27 

1.40±0.18 

- 470±100 

-i-430±190 

Ovilf 

22.101 

8.96±0.55 

1.29±0.24 

- 250±20 

-i-360±160 

N VII Lya 

24.779 

3.69±0.94 

0.80±0.16 

-H 15±360 


C VI Lya 

33.734 

4.69±0.68 

1.21±0.31 

- 300±120 

- 860±330 

Cvf 

41.472 

6.44±1.99 

1.01±0.40 

- 380±160 



4.2. Line Shifts 

We used the same set of lines to look for redshifts or blueshifts. 
First of all we find for most lines of the central region a dif¬ 
ference between the measured positive and negative spectral 
order position of about 0.02 A. One would however expect a 
difference of ^ 0.01 A, based on the difference between the 
peak position of the Gaussian fit to the zero order distribution 
and the zero point of the wavelength scale. The underlying as¬ 
sumption is that the line radiation distribution coincides with 
the total zero order radiation distribution. We verified the line 
wavelength obtained above by fitting the measured line shape 
to a Gaussian and forced the line width to be equal to the zero 
order width; that has very little effect on the line wavelength 
determination. The results of fitting the positive and negative 
spectral order lines together, are summarized in Table Q. We 
then looked for possible differences in redshift or blueshift be¬ 
tween the Secondary and Primary region. 

Fitting for the Secondary region, the positive and negative 
spectral orders together lead to redshifts, also listed in Table 
(All shifts are calculated with respect to the systemic velocity). 

5. Long Wavelength L-shell Line Detections 


Given the limited and not homogeneously tested accuracy 
of the wavelengths, and also the possibility of some blending 
we cannot always be sure of the derived fluxes in each indi¬ 
vidual case. However, we think that the set of fluxes still re¬ 
flects the overall behaviour of these ions. Not all the fluxes 
are statistically significant, but for the weakest lines the value 
given plus its nominal error bar gives at least an order of ma^ 
nitude estimate of the column density of the parent ion. Table g 
lists the ionization parameter, for which the ion concentra¬ 
tion reaches a maximum. We have also made model predic¬ 
tions for the line fluxes. We assume that for each line the emit¬ 
ted flux is equal to the number of photons taken away from the 
power law continuum (F = 1.7, covering factor times luminos¬ 
ity fLx = 10"^® erg/s, cf. the RGS analysis of Kinkhabwala et 
al. 2002a). To be more specific, we approximate the equivalent 
width for the absorption line by; 


W = 


AtTv 


OO 

y [1 - exp(-Toe“2^ 


( 1 ) 


with To equal to the optical depth of the line at the line center, 
given by 


We have also analysed the long wavelength part of the LETGS 
spectrum of the Primary component (beyond 30 A). This part of 
the spectrum contains the L-shell complexes of Ne, Mg, Si, S, 
Ar and Ca. We have determined the fluxes of the strongest res¬ 
onance lines in the 15 — 60 A band (see Table ^). Wavelengths 
and oscillator strengths for these lines are the same as those 
quoted in Kaastra et al. (2002). The data are taken from what 
is used in their slab and xabs models, and for the relevant 
ions they are derived mostly from the compilation by Verner 
et al. (1996); updates in some cases, using calculations with 
HULLAC (Bar-Shalom et al. 2001), were provided by Liedahl 
(see Kaastra et al. 2002 for more details). In addition, we have 
included a number of Fe-L transitions in the table. For some 
ions we did not use the strongest line. An example is Fe XlX, 
where the strongest lines are blended by lines from NelX. 
Finally, in a few cases we co-added lines from the same ion 
if they could not be resolved with the LETGS. 


TO = O.lOe/A^ioA/avpoo. (2) 

Here / is the oscillator strength, A the wavelength in A, Uv.ioo 
the velocity dispersion in units of 100 km/s and the column 
density of the ion in units of 10^® cm“^. 

In the above, we neglected the damping wings of the lines. 
Equation [| gives the equivalent width for an absorption line. 
The power law continuum is given by dN/dX = 8.67 x 
10"^® A“° ® with A in A. The predicted emission line flux is then 
simply F = WdN/dX. We have divided the 40 emission lines 
into 4 groups of ten, each spanning a range in ionization pa¬ 
rameter. Eor each group, we have determined the velocity dis¬ 
persion (Tv and total hydrogen column density A^h that best de¬ 
scribe the observed line fluxes (assuming Solar abundances). 
The results are summarized in Table 

All hydrogen column densities are close to the value of 
logA^H = 21.9, and we have adopted that value uniformly 
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Table 5. Wavelengths, oscillator strength, ionization parame¬ 
ters and line fluxes for the strongest lines of each species in 
the Primary component. Line strengths are given in units of 
10^® ph/s. 


Ion 

A 

(A) 

/ 

logC 

Flux 

observed 

Flux 

model 

NeVIII 

57.75 

0.01 

1.0 

2.7± 4.2 

2.2 

Mg IX 

48.34 

0.14 

1.2 

3.6± 1.3 

3.1 

MgX 

57.88 

0.22 

1.5 

6.0± 5.0 

4.7 

Si VIII 

50.02 

0.12 

0.6 

2.9± 1.4 

2.9 

Si IX 

55.30 

0.95 

0.9 

5.9± 3.0 

4.5 

Six 

50.52 

0.60 

1.3 

4.3± 1.7 

5.1 

Si XI 

43.76 

0.48 

1.7 

4.7± 1.0 

3.2 

Si XII 

40.91 

0.23 

2.0 

3.0±5.1 

2.2 

Svil 

51.81 

0.61 

0.0 

3.4± 1.8 

3.9 

S VIII 

52.77 

0.89 

0.4 

3.7± 2.3 

4.1 

Six 

47.43 

0.38 

0.8 

2.7± 1.0 

3.2 

Sx 

42.54 

0.67 

1.1 

4.0± 4.3 

3.3 

Sxi 

39.24 

1.05 

1.4 

3.3± 2.2 

4.1 

Sxii 

36.40 

0.62 

1.7 

3.1± 1.3 

2.5 

Sxiii 

32.24 

0.45 

2.0 

2.5± 0.8 

2.1 

Sxiv 

30.44 

0.34 

2.4 

1.6± 1.0 

0.6 

Arx 

37.43 

0.65 

1.0 

3.6± 1.7 

1.4 

ArXI 

34.33 

0.42 

1.3 

2.9± 1.0 

1.0 

ArXII 

31.39 

0.87 

1.5 

4.2± 0.9 

1.5 

ArXIII 

29.32 

0.91 

1.8 

1.7± 0.7 

1.2 

ArXIV 

27.47 

0.65 

2.1 

0.7± 0.5 

0.8 

CaXI 

30.45 

2.34 

1.2 

1.6± 1.0 

1.5 

CaXII 

27.97 

1.30 

1.3 

0.2± 0.4 

1.0 

CaXIII 

25.53 

0.47 

1.6 

0.3± 0.2 

0.4 

CaXIV 

24.13 

0.89 

1.9 

0.4± 0.2 

0.6 

Caxv 

22.73 

0.93 

2.1 

0.5± 0.2 

0.4 

CaXVI 

21.45 

0.66 

2.4 

0.9± 0.5 

0.2 

CaXVII 

19.56 

0.65 

2.5 

0.2± 0.1 

0.2 

CaXVIII 

18.70 

0.37 

2.8 

0.3± 0.3 

0.1 

Fe XIV 

58.96 

0.28 

1.4 

5.4± 8.5 

4.5 

Fe XV 

52.91 

0.41 

1.6 

4.3± 2.5 

4.7 

Fe XVI 

50.35 

0.16 

1.6 

5.4± 1.9 

2.4 

Fe XVII 

15.01 

2.72 

2.1 

1.6± 0.3 

1.6 

Fe XVIII 

14.20 

0.88 

2.3 

0.8± 0.2 

1.2 

Fe XIX 

13.80 

0.21 

2.5 

0.4± 0.2 

0.3 

Fe XX 

12.84 

0.93 

2.8 

0.7± 0.2 

0.4 

FeXXI 

12.29 

1.24 

3.0 

0.3± 0.1 

0.4 

FeXXII 

11.71 

0.67 

3.1 

0.2± 0.2 

0.3 

Fe XXIII 

8.30 

0.14 

3.3 

O.liO.l 

0.1 

Fe XXIV 

10.64 

0.37 

3.4 

0.2± 0.2 

0.2 


Table 6. Velocity dispersion and column densities for different 
ranges in log 


log^ 

CTv 

range 

log Ah 

range 

0.0-1.2 

170 

120-220 

21.9 

lower limit 

1.3-1.6 

220 

130-300 

21.9 

21.8-22.1 

1.7-2.3 

150 

120-210 

21.8 

21.5-22.0 

2.4-3.4 

40 

30-60 

22.0 

21.5-22.7 


in the predicted line fluxes of Table ^ For the velocity dis¬ 
persions in that table we have adopted the values of Table ^ 
instead. We find evidence for a significantly lower velocity dis¬ 
persion at higher values of the ionization parameter. This may 


explain partly the problem we had with the H-like Ne and Mg 
lines, where we had fixed CTv to 100 km/s. The tendency of a 
smaller velocity dispersion for the most highly ionized mate¬ 
rial is also found in the Seyfert 1 galaxy NGC 5548, where the 
high ionization line of Si XIV Lya is consistent with a velocity 
component that in the corresponding UV spectra shows up as 
much narrower than the lower ionized lines. 


6. Summary and Conclusions 


We have shown that the Chandra LETGS spectrum of 
the Primary region (comprising the bulk of the flux from 
NGC 1068) is consistent with emission from a photoionized 
and photoexcited cone of warm plasma, which is in accord 
with previous observatio ns of th is object with XMM-Newton 
RGS (Kinkhabwala et al. 2002b). We have also shown, using 
the spatial resolution capability for angles of a few arcsec of 
the Chandra grating spectrometers, that the X-ray spectrum 
of a separate region (4" to the NE), which we refer to as the 
Secondary, can be explained by a similar model, except with 
a factor of a few lower inferred radial ionic column densities. 
Any additional emission in either spectrum due to hot, colli- 
sional plasma is negligible. 

The lower radial ionic column densities of the Secondary 
region as compared with the Primary region are naturally ex¬ 
plained in the context of an irradiated cone of plasma. The 
radial ionic column density through each region is given in 
Kinkhabwala et al. (2002b) as; 


~ 1.6 X 10^4Ap 


9 


frn 


r /.] 

[ (/Ax) 1 

r e 1 

0.5 

1043 gj.gg g-1 

1 ergs cm 


cm 


( 3 ) 


where Az is the abundance relative to hydrogen, g is the radial 
filling factor, / = 12/471 is the usual covering factor, rmin.pc is 
the minimum radius of the cone from the nucleus in pc, and fi 
is the fractional ionic abundance. The radial dependence of r~^ 
in this equation provides a simple explanation for the observed 
factor of a few drop in column densities in the Secondary region 
versus those inferred for the Primary region (Table H), which is 
closer to the nucleus by at least a similar factor of a few. 

An estimate of the radial filling factor of the plasma in the 
Secondary region is also made possible by our knowledge of 
rmin (~4"; 288 pc) for this region. The inferred radial column 
density of iV/Jj/ = 2.10^^ for Ovil in the Secondary region, 
along with Az = 4.89 x 10“^ and / < 0.1, imply a radial 
filling factor of p < 0.07 for plasma at the typical ionization 
parameter for O VII of ^ = 10 ergs cm s“^. 

The remarkable overlap of the [O III] and X-ray cones ap¬ 
pears to indicate that these regions are related (Young et al. 
2001). X-ray spectroscopy of the Primary region reveals excess 
blueshifts, whereas spectroscopy of the Secondary region re¬ 
veals excess redshifts. These blueshifts and redshifts are com¬ 
parable to the shifts observe d in the same respective regions for 
optica l/UV lines (Cecil et al. 199C , 2002 ; Crenshaw & Kraemer 
|2000 ), kinematically linking the X-ray emission regions with 
lower ionization state emission regions. In addition, Kraemer 
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& Crenshaw (2000) argue that the optical/UV emission in the 
cone is also ultimately powered by photoionization due to the 
inferred nuclear continuum. 

Using both K-shell and L-shell line series, we find that 
the relative abundances of all observed elements are consistent 
with Solar abundances, aside from N. From the relative strength 
of the column densities derived from K-shell emission features, 
we hnd evidence for a relative factor of 2—3 more N compared 
with C or O. Further investigation of elemental abundances is 
left for the future. 

The K-shell transitions as well as the L-shell transitions for 
multiple ions, imply that a broad range in ionization parameter 
is present. The long wavelength emission line spectrum is very 
similar to what is found in Seyfert 1 galaxies in absorption. 
These observations provide further conhrmation of the unihca- 
tion schemes, wherein Seyfert 1 and Seyfert 2 galaxies consti¬ 
tute the same objects seen from different viewing angles. The 
data also suggest a variation of velocity dispersion as a function 
of ionization parameter (see Table |^. 
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